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BENTHIC COMMUNITY STRUCTURE IN TWO ADJACENT 
STREAMS IN YELLOWSTONE NATIONAL PARK FIVE YEARS 
AFTER THE 1988 WILDFIRES 


G. Wayne Minshall!, Christopher T. Robinson!, Todd V. Royer!, and Samuel R. Rushforth? 


ABSTRACT.—Physical characteristics, benthic macroinvertebrates, and periphyton assemblages in two adjacent head- 
water streams in Yellowstone National Park were evaluated five years after the 1988 wildfires. The catchment of one 
stream was burned by wildfire (burned stream) while the other catchment was unburned (unburned stream). Physical 
measures revealed channel alteration in the burned stream relative to the unburned stream. Periphyton biomass was 
lower in the burned than the unburned stream (29.2 vs. 50.5 gin? AFDM, respectively), further demonstrating the 
unstable physical conditions of that system. Kendall's coefficient of concordance (an index of similarity) between diatom 
assemblages was 0.22, indicating distinct assemblage composition between streams. Navicula permitis Hust. was the 
most abundant diatom in the burned stream while Hannaea arcus (Ehr.) Patr. was dominant in the unburned stream. 
Macroinvertebrate taxa richness, density, and biomass were all greater in the unburned stream, although Chironomidae 
was the most abundant taxon in both streams. Results suggest the removal of terrestrial/riparian vegetation by wildfire 
can directly influence stream benthic assemblages by altering the inherent disturbance regime of the physical habitat 
templet. 
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Physical disturbance, acting at various spa- 
tial and temporal scales, often is the predomi- 
nant factor structuring stream benthic com- 
munities (Minshall 1988, Resh et al. 1988). 
Further, physical disturbances may be viewed 
in a hierarchical framework, with the effects of 
small-scale disturbances altered (intensified or 
mediated) by large-scale disturbance events 
(sensu O Neill et al. 1986). Wildfire, as a large- 
scale disturbance, directly influences stream 
biotic structure and function by affecting the 
physical habitat of stream ecosystems (Minshall 
et al. 1989, Minshall and Brock 1991, Richards 


and Minshall 1992, Robinson et al. 1994, Mihuc 
et al. in press, Robinson and Minshall in press). 
In lotic ecosystems, physical disturbance also 
may constrain the establishment of biotic con- 
trols, such as competition and predation, on 
benthic community structure (McAuliffe 1984, 
Minshall and Petersen 1985, Resh et al. 1988). 

Wildfire burned extensive portions of the 
Greater Yellowstone Ecosystem during the 
summer of 1988. Over 32% of the streams in 
Yellowstone National Park (YNP) were affected 
to varying degrees by wildfires (Minshall et al. 
1989, Minshall and Brock 1991, Robinson et al. 
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199-4). Minshall and Brock (1991) summarized 
the immediate effects of the fires on YNP's 
stream ecosystems and hypothesized on the 
mid-term (10-25 yr) and long-term (50-300 vr) 
elfeets. They suggested that most adverse 
short-term cffeets on streams resulted from 
increased sediment load and channel erosion 
caused by increased overland runoll following 
precipitation events and snowmelt. The in- 
tensity and frequency of short-term effects 
were hypothesized to decrease by vear 5 as 
riparian conditions improve (see Richards and 
Minshall 1992). In general, mid- and long- 
term effects on streams, including reeovery to 
prefire conditions, should correspond to vegeta- 
tive regrowth in burned catchments (Minshall 
et al. 1989, Minshall and Brock 1991). 

The present study compared physical char- 
acteristics and benthic community structure in 
two streams five vears after the 1988 wildfires. 
The streams are adjacent second-order (after 
Strahler 1952) tributaries of the South Fork 
Cache Creek. The catchment of one stream 
was burned during the 1988 wildfires, while the 
catchment of the other was essentially un- 
burned. The spatial arrangement of these 
streams (adjacent basins) provided a treatment/ 
reference situation where confounding factors 
of climate and geology are minimized when 
comparing differenees among the study 
streams. However, the study lacks true repli- 
cation of the burned and unburned treatments 
(sensu Hurlbert 1984) and must be viewed as 
a simple comparison study. Phenomenologieal 
studies and/or two stream comparisons are 
common in stream ecology (e.g., Wallace et al. 
1986, Robinson et al. 1993, Scarsbrook and 
Townsend 1993) and are eapable of providing 
valuable insights (Townsend 1989). The present 
study can be viewed as a natural “experiment” 
with observed differenees between the two 
streams attributed to the effects of wildfire. In 
that context, the study provides insights on 
general patterns of lotic ecosystem recovery to 
an unpredietable, large-seale disturbance 
(Townsend 1989, Lamberti et al. 1991). 


METHODS 


The study streams, located in the northeast 
corner of YNP. were surveyed on 19 July 1993. 
One stream had over 80% of its catchment 
burned during the 1988 Yellowstone wildfire 
(hereafter, burned stream; 110°01'30’W, 
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44°50'00"N), while less than 10% of the eateh- 
ment of the other stream was burned (hereafter, 
unburned stream; 110°01’00"W, 44°49'30"N). 
Climate of the area is typical of the northern 
Roeky Mountains, with precipitation primarily 
oecurring as snow during the winter months. 
Both streams drain eatchments primarily veg- 
etated (prior to the fire in the burned stream) 
by coniferous forests of lodgepole pine (Pinus 
contorta) and Engelmann spruce (Picea engel- 
mannii). Riparian vegetation consisted of wil- 
low (Salix), rose (Rosa), and alder (Alnus). 

Surveys were conducted approximately 0.5 
km above the confluence of the two streams. 
Physieal characteristics were measured in each 
stream at five cross-sectional transeets, each sit- 
uated approximately 50 m apart. Measurements 
made at each transeet included stream width 
at baseflow, stream width at bankfull discharge, 
and stream eross-sectional profile (for caleula- 
tion of width:depth ratios). Discharge was eal- 
culated in each stream at the most suitable 
transect following the methods of Platts et al. 
(1983). In addition to measurements at each 
transect, 100 randomly seleeted rock substrata 
along a 100-m length of stream (loeated within 
the outermost cross-sectional transeets) were 
measured for size (length of the longest axis) 
and pereent embeddedness. Embeddedness 
was defined as the pereent eoverage of the 
rock (three-dimensional surface) by fine sedi- 
ments. Large boulders that protruded through 
the water surface were not used in substratum 
size measurements. Water depth and near-bed 
water velocity also were reeorded at each of 
the 100 random locations. Near-bed water 
velocity was measured with a small Ott C-1 
eurrent meter approximately 2 cm above each 
substratum, 

One periphyton sample was eolleeted from 
a suitable (flat-surfaced, medium-sized) rock 
substratum at each cross-sectional transeet 
using a method described in Robinson and 
Minshall (1986). Samples were frozen in the 
field in a Taylor-Wharton 3DS dry shipper 
eharged with liquid nitrogen and returned to 
the Jaboratory for processing. In the laboratory, 
samples were extracted in 10 ml of methanol 
for 24 h (Holm-Hansen and Riemann 1978). 
One 3-ml subsample was then removed from 
each sample and analyzed for chlorophyll a 
using a Gilford Instruments (Model 2600) spec- 
trophotometer. The remaining periphyton 
material from each sample was used for algal 
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biomass determination, expressed as grams 
ash-free dry mass (AFDM) per m?. The mater- 
ial was dried at 50°C for 24 h, weighed on a 
Sauter balance (Model AR 1014), ashed at 
550°C for a minimum of 3 h, rehydrated, 
redried at 50°C, then cooled to ambient tem- 
perature in a desiccator and reweighed. The 
difference in weights equaled the AFDM of 
the sample. 

Diatom samples were collected in each 
stream, after Robinson and Rushforth (1987), 
from three to five rock substrata representing 
the predominant habitat type (typically riffles). 
Samples were composited, preserved with 5% 
formalin, and returned to the laboratory. The 
composite sample was boiled in concentrated 
nitric acid, rinsed, mounted in Naphrax moun- 
tant, and examined under 1000X oil immersion 
using a Zeis RA microscope with Nomarski 
optics (St. Clair and Rushforth 1976). Relative 
abundances of diatom taxa were determined 
by counting a minimum of 1000 diatom valves 
from each stream. Diatoms were analyzed in 
terms of species richness, Simpson's index, 
and Kendall's coefficient of concordance (an 
index of similarity using all taxa with a relative 
abundance >1%). Other algal groups such as 
Chlorophyta (green algae) and Cyanobacteria 
(blue-green algae) were not abundant at the 
time of sampling and thus were not consid- 
ered in the present study. 

One benthic sample was collected from a 
riffle/run habitat (pools were rare and not sam- 
pled) near each transect and analyzed for 
macroinvertebrates and benthic organic mat- 
ter (BOM). Samples were collected using a 
Surber sampler (250 um mesh), preserved 
with 5% formalin, and returned to the labora- 
tory. Woody debris >5 cm in length that was 
collected in the benthic samples was rinsed of 
invertebrates and removed from the samples. 
In the laboratory, macroinvertebrates were 
hand-sorted from the benthic detritus with the 
aid of a 3X dissecting microscope, identified to 
the lowest feasible level (usually genus), enu- 
merated, dried at 50°C for a minimum of 48 h, 
then cooled to ambient temperature in a des- 
iccator and weighed. Dry weights, in milli- 
grams, were determined on a Cahn (Model 
25) electrobalance. The benthic detritus from 
each sample was used for BOM determination. 
The quantity of BOM, expressed as g AFDM/ 
m2, was determined as described above for 
periphyton. Macroinvertebrates were analyzed 
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in terms of density (no./m?), biomass (mg/m2), 
taxa richness, Simpson's index, and relative 
abundances. 

Chi-squared analysis was used to test for 
statistical differences in median substratum 
size between the two streams (Zar 1984). 
Independent sample t tests were used to com- 
pare the other characteristics for differences 
between the two streams. Prior to the t-test 
analysis all data were log (x + 1) transformed, 
except substratum embeddedness and the rela- 
tive abundance of invertebrate taxa (both per- 
centage measures), which were arcsine (square 
root [x]) transformed (Zar 1984). Tabular re- 
sults are presented as untransformed means 
and standard deviations. All statistical analyses 
were performed on SYSTAT (Wilkinson 1990). 


RESULTS 


Baseflow discharge was equal in the two 
streams (0.2 m/s), reflecting the similar catch- 
ment size of the bumed (22 km?) and unbumed 
(26 km?) streams. Mean baseflow width, near- 
bed water velocity, and BOM were not signifi- 
cantly different between the two streams (P > 
.05). Substratum embeddedness was signifi- 
cantly greater in the burned stream (P = .01), 
although the difference between mean values 
was not large (burned = 62.9, unburned = 
52.8). It is not known whether this statistical 
difference was biologically meaningful or sim- 
ply a reflection of the large sample size (n = 
100). 

Water depth at baseflow (P < .01) was lower 
and stream width at bankfull discharge greater 
(P = .03) in the bummed stream than the un- 
burned stream. Although not statistically sig- 
nificant (P = .06), the ratio of stream width: 
depth was greater in the burned than the un- 
burned stream (216 and 91, respectively). The 
general appearance of the two streams was dis- 
tinctly different (Fig. 1; Minshall personal ob- 
servation). Large, woody debris and streamside 
riparian vegetation, which provide bank and 
channel stability, were noticeably absent in the 
burned stream. 

Mean substratum size was not significantly 
different between the two streams (P > .05) in 
1993, possibly because large boulders were 
not recorded in the measurements (see Fig. 1). 
We collected additional data on substratum size 
in August 1994 and included large boulders in 
the measurements. Further, substrata within 
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Fig. 1. Representative photographs of the burned (upper) and unburned (lower) streams five years after the 1988 
wildfire. Note absence of large, woody debris and streamside riparian vegetation in the burned stream. 
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the bankfull channel were measured in 199-4, 
in contrast to measures being reeorded only 
within the baseflow channel in 1993. The 199-4 
results showed that mean substratum size was 
significantly larger in the unburned than in 
the burned stream (P < .01). 

A comparison of median substratum size 
showed similar results to that of mean substra- 
tum size. Median substratum size was not dif- 
ferent between the two streams when mea- 
surements excluded large boulders and were 
confined to the baseflow ehannel (P > .05). 
However, when measurements inchided large 
boulders and encompassed the bankfull chan- 
nel, the difference in median size was signifi- 
eant (P < .01). Whether large boulders were 
present in the burned stream prior to the wild- 
fire has yet to be determined. However, in other 
streams influeneed by intensive wildfire, large 
boulders were observed to be buried by inor- 
ganie debris (primarily gravel and fine sedi- 
ments) within five years following wildfire 
(Minshall personal observation). 

The burned stream contained less periphy- 
ton chlorophyll a (P = .06) and AFDM (P < 
01) than did the unburned stream (Table 1). 
Diatom species richness was greater in the 
burned (34 taxa) than in the unburned stream 
(27 taxa; Table 2). Simpson's index was lower for 
the burned than the unburned stream (0.12 
and 0.42, respectively). Kendall's coefficient of 
concordance for the two diatom communities 
was 0.22, suggesting distinet assemblage com- 
position among sites. For example, Navicula 
permitis Hust. was the most abundant species 
in the burned stream, constituting 24.7% of 
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the assemblage, while Hannaea arcus (hr) 
Patr. comprised 63.1% of the assemblage in 
the unburned stream (Table 2). 

Mean macroinvertebrate density and bio- 
mass were lower in the burned than unburned 
stream (Table 3), but the differences were not 
significant (P > .05). For example, mean den- 
sity in the burned stream was 9960 individuals/ 
m2, while the unburned stream had 16,950 indi- 
viduals/m?, and mean biomass (dry weight) 
was 1960 and 3200 mg/m? in the burned and 
unburned streams, respectively. Taxa richness 
and Simpsons index both were reduced in the 
burned stream, although the differenee was 
significant only for Simpson's index (P = .0-4) 
(Table 3). The burned stream contained a 
mean of 15 taxa per benthic sample compared 
to a mean of 20 taxa for the unburned stream. 
The mean Simpson's index was 0.57 for the 
burned stream and 0.73 for the unburned 
stream. Chironomidae was the most abundant 
taxon in both streams (Table 4), although their 
relative abundance was significantly greater (P 
= .03) in the unburned stream. There were no 
statistical differenees (P > .05) in relative 
abundances of other taxa common to both 
streams (Hydraearina, Simuliidae, Baetis bi- 
caudatus, Cinygmila, and Zapada coliunbiana). 


DISCUSSION 


Alterations of the surrounding terrestrial 
landseape by major unpredictable disturbances 
such as hurricanes, voleanie eruptions, or wild- 
fire directly influence streams draining the 


TABLE 1. Means (SD) and P values for physical characteristics measured in the study streams. 








Burned Unburned P value 
Baseflow width (m) 5.9 (2.6) to (1S) 46 
Near-bed velocity (cim/s) 10.3 (0.1) Oz 40:1) rf 
BOM (g/m2) 15 (0.8) 25 Ie?) .20 
Embeddedness (%) 62.9 (28.5) 52.0 (309) 01 
Baseflow depth (cm) 16:9.2e0( 11.0) OA le) = 00 
Bankfull width (m) 35.0 (11.9) 16.6 (1.4) 03 
Bankfull width:depth ratio 216 = (101) 9) Pal) .06 
Periphyton chl @ (mg/m2) &:9) (5.5) 82.1 (19'5) .06 
Periphyton AFDM (g/m?) 29:2 (3.7) 50.5 (9.3) <.01 
“Mean substratum size (cm) ee 410.5) 16.0; (14-2) 26 
Mean substratum size (cm) leyome lo.) 27.6 (27.4) <.0] 
Median substratum size (cm) 12.0 11.0 fie! 
Median substratum size (cm) 11.0 17.0 <.01 





: Baseflow channel, large boulders excluded. 
Bankfull channel, large boulders included. 
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fABLE 2. Community measures and relative abundances 
C) for the diatom assemblage of each study stream. 





Burned Unburned 

Se rN ee ee 
Species richness 34 27 
Simpson's index (C) 0.12 0.42 
Navicula permitis lust. 24.7 Dal 
Nitzschia disstpata 

(Kuetz.) Grun. 17.3 4.2 
Achnanthes lanceolata 

(Breb.) Grun. 9.8 Ley 
Nitzschia paleacea Grun. 7.6 3. 
Navicula arvensis Hust. 4.5 2.4 
Hannaca arcus (Ehr.) Patr. 2A 63.1 





affected watersheds. For example, the Mt. St. 
Helens eruption of 1980 dramatically changed 
drainage patterns and river networks, elimi- 
nated terrestrial vegetation, and caused major 
debris flows that scoured stream channels 
(Wilzbach et al. 1983, Hawkins 1988). However, 
high spatial variation in the intensity of these 
major disturbances may occur, causing tempo- 
ral differences in recovery patterns (Yount and 
Niemi 1990). In catchments of YNP the rela- 
tive area burned ranged from <10% to >90% 
(Minshall and Brock 1991). Further, the 
degree of alteration of stream habitat was 
highly correlated with percent of catchment 
burned (Robinson and Minshall in press). 

In the present study significant differences 
were observed in the benthic habitat of the two 
streams. The width:depth ratio of the burmed 
stream was greater than that of the unburned 
stream. Anderson (1992) also observed in- 
creased width:depth ratios following major 
disturbances in streams of the Cascade Moun- 
tains. With large boulders included in the mea- 
surements, the unburned stream exhibited sig- 
nificantly greater substratum size. Gurtz and 
Wallace (1984) demonstrated that Jarge sub- 
strata could mediate the effects of large-scale 
disturbances by providing stable habitat for 
benthic organisms. At the time of sampling, 
the burned stream did not contain the larger- 
sized substrata found in the unburned stream. 
It is probable that the larger substrata in the 
burned stream were buried by inorganic sedi- 
ments following the wildfire, as has been 
observed in other YNP streams (Minshall per- 
sonal observation). Thus, one effect of the wild- 
fire appeared to be alteration of the substrata 
in such a manner as to make the benthic habitat 
more susceptible to future disturbances (e.g,, 
Gurtz and Wallace 1984). 
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TABLE 3. Means (SD) and P values of macroinvertebrate 
community measures for the study streams. 


Burned Unburned P value 
Density (no./m?) 9963 (4730) 16,945 (8995) al 
Biomass (mg/m?) {956 (1056) 3198 (1274) 19 
Taxa richness 150537) 20.4 (2.6) 58 


Simpson's index (C) 0.57 (0.09) 0.73 (0.10) 04 


Lamberti ct al. (1991) found that faunal 
densities and macroinvertebrate species rich- 
ness had recovered within one year following 
a major debris flow in an Oregon stream. In 
central Idaho, however, streams disturbed by 
wildfire and unburned reference streams 
showed little similarity in macroinvertebrate 
assemblages, even after five years of recovery 
(Richards and Minshall 1992). Similarly, in the 
present study the influence of wildfire was still 
apparent after five years. Macroinvertebrate 
community structure was not similar between 
the two systems, despite their close proximity 
to each other (0.5 km). Most researchers agree 
that recovery of the benthic community will 
correspond to recovery of the surrounding 
landscape (Steinman and Lamberti 1988, 
Minshall et al. 1989, Lamberti et al. 1991, 
Minshall and Brock 1991, Anderson 1992, 
Richards and Minshall 1992, but see Hawkins 
1988). 

Primary producers (lotic algae) may recover 
sooner than consumers (macroinvertebrates and 
fish) because of their much shorter life cycles, 
and subsequently they may influence recovery 
of the higher trophic levels (Steinman and 
MeIntire 1990). In the present study, peri- 
phyton biomass (as AFDM) in the unburned 
stream was I.7X greater than in the burned 
stream, implying a present lack of recovery by 
primary producers in the burned system. 
Macroinvertebrate taxa richness also was 
greater in the unburned stream than in the 
burned stream. How functional or structural 
recovery of macroinvertebrates is related to 
algal recovery following wildfire has yet to be 
determined, but provides an interesting and 
important avenue for future research. Algae 
have shorter life eyeles and reduced mobility 
relative to macroinvertebrates, and possibly 
the two groups respond differently to large- 
scale disturbances. 
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TABLE 4. Mean (SD) density and relative abundance of 
the 10 most abundant macroinvertebrate taxa from each 
stream. These taxa constituted >90% of the assemblage in 
their respective streams. 








Taxa Density Relative 
(no./m2) abundance (%) 
222 ..------------- Burned - ----------------- 
Chironomidae 5437 (2218) 59.7 (19.2) 
Simuliidae 27371 (583) 18.2 (24.1) 
Baetis bicaudatus 576 (3-46) 8.1 (6.3) 
Hydracarina 148 (232) We (2.3) 
Epeorus albertae 95 (123) LEO al 2 
Zapada columbiana 90 (72) 1.0 (0.7) 
Epeorus longimanus SO (125) 0.5 (0.7) 
Cinygmula fo 7) 0.8 (0.7) 
Rhithrogena Si\p2) 05 x03) 
Ameletus cooki 22S) U3 (0.2) 
~---------------- Unburned - - --------------- 
Chironomidae 14,676 (S259) 84.7 (6.7) 
Hydracarina 372 (253) 2.) (1:6) 
Cinygmula 314 (187) 1.7 (0.8) 
Zapada columbiana 310 (215) 1.6 (0.6) 
Simuliidae 307 (181) 2.5 (5.3) 
Baetis bicandatus 125 (60) ee 7) 
Drunella coloradensis 75 (64) 0.4 (0.5) 
Rhyacophila angelita >) 0.4 (0.2) 
Rhyacophila tucula 5 (29) 0.3 (0.1) 
Kogotus 520(35) 0.3 (0.2) 


After five years of recovery, the channel of 
the burned stream still appeared unstable as 
indicated by different diatom assemblages 
between the two streams. For example, the 
small, adnate diatom Navicula permitis Hust. 
was predominant in the burned stream but 
was found in relatively low abundance in the 
unburned stream. N. permitis was predomi- 
nant in other YNP streams influenced by the 
1988 wildfires, and it has been suggested that 
a diatom community with an abundance of N. 
permitis is indicative of more physically dis- 
turbed stream environments (Robinson et al. 
1994). Further, Robinson et al. (1994) showed 
diatom recovery among, 14 streams in Yellow- 
stone was inversely related to degree of dis- 
turbance by wildfire. Similarly, Steinman and 
Lamberti (1988) found little recovery, after six 
vears, in the composition of algal communities 
in intensively disturbed streams of Mt. St. 
Helens. In summary, benthic community re- 
covery patterns appeared to be related to the 
recovery of stream physical habitat which, five 
years after the 1988 wildfires, still displayed 
evidence of instability. When examining the 
recovery of benthic communities following 
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large-scale disturbance, one must remain aware 
of the connections between the terrestrial land- 
scape, lotic habitat, and benthic organisms. 
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